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Abstract 



^ ■ Non-positive definite, global inviscid invariants similar to helicity are 

b : 

discussed for two types of shell models and evidence for a new role for 
helicity in Navier-Stokes turbulence is presented. It is suggested that 
the extra invariants play the role of triggering the intermittent cascade 
of energy to small scales characterized by pulses. These invariants also 
determine where a transition to chaos appears. New analysis of numer- 
ical experiments with existing models is suggested and a new class of 
shell-models where the dynamical interactions of a second quadratic in- 
variant are closer to those of helicity in the Navier-Stokes equations is 
introduced. The place of the popular GOY model within this class is 
discussed. PACS number 47.27. Eq 
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1 Introduction 



One of the outstanding questions in our understanding of fully-developed tur- 
bulence is the mechanisms by which the cascade of energy to small scales is 
maintained. That the cascade is intermittent is well recognized, but the phe- 
nomenological and dynamical models used to address the problem are rarely 
connected to the dynamics in the full Navier- Stokes equations. In this note 
it is shown that a popular model for explaining turbulent intermittency, the 
GOY model fl], g], shares some symmetries with terms in a decomposition of 
the spectral Navier-Stokes equations into the interactions between its helical 
components ||. The essential common property that is identified in both the 
GOY model and in Navier-Stokes is the importance of interactions between 
components with oppositely signed helicity. This new role for helicity sup- 
porting an intermittent cascade contrasts strongly with helicity's previously 
identified role in blocking the cascade j|, ||, a role that was proposed earlier 
based on an analogy to the way magnetic helicity creates force- free states [||. 

In the GOY model with the standard parameters for three-dimensional tur- 
bulence, interactions between modes with oppositely signed helicity occur nat- 
urally as the sign of helicity reverses between neighboring shells |7| . Whether 
Navier-Stokes turbulence follows a similar path is a more difficult question be- 
cause there are several paths, characterized by interactions between different 
components of the helicity, that the cascade can follow. One way to investigate 
this question is to consider several variants of the GOY models that investigate 
each path in the full Navier-Stokes equations individually or in unison. One 
question would be how strongly the statistical behavior of the cascade depends 
on the symmetries in the different models. Another line of investigation is to 
determine which path the cascade follows in the full Navier-Stokes equations. 
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In this note, preliminary results following both of these approaches are pre- 
sented and the direction of a more complete study is presented. Also included 
will be new analysis of the Kerr-Siggia shell model ||, which has a cubic, non- 
positive definite invariant that shares some dynamical properties with helicity 
and served as an inspiration for part of this new proposal. In all of these 
cases it will be argued that competition between the transfer of energy and the 
transfer of the generalized helicity could explain the presence of numerically 
observed chaotic dynamics and intermittency in the energy cascade. For the 
GOY models, new results on intermittency corrections will be used to illustrate 
the importance of the second quadratic invariant in the energy cascade. 

In section 2, the new results for the Kerr-Siggia model will be reported. In 
section 3, the GOY model is reviewed and how the inviscid conserved quan- 
tities and the dynamics depend on the free parameters present in the model 
is discussed. In section 4, an argument that predicts the transition from a 
trivial dynamics (dominated by the presence of an attractive fixed point) to a 
fully chaotic regime for some critical values of the free parameters is discussed. 
Some numerical results for the energy transfer are also discussed. In section 
5, new versions of the GOY model are introduced by considering explicitly the 
possibility of having shells which transport positive or negative helicity exactly 
as occurs in the Navier-Stokes equations. Two preliminary calculations with 
the full Navier-Stokes equations that support the importance of the interac- 
tions between components with oppositely-signed helicity are presented. Some 
problems that can be studied by using the new variant of the GOY model are 
discussed and some new analysis of the full Navier-Stokes equations that could 
be done to illuminate these properties is presented. 
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2 Pulse scaling of Kerr-Siggia 

The Kerr-Siggia model || is a shell model with one complex variable per shell, 
originating from a decimation of the possible interactions between triads in 
Burgers equation. With a simple forcing and eddy viscosity the equations were 

dui/dt = e/u[ + 1iuiu\ 

du n /dt = kn-iiiu 2 ^ + 2u n+1 u* n ) (1) 

du^/dt = kN-ii(u 2 N _ 1 + v e \uN\uN) 

where e was the average energy input and dissipation, v e was taken to be 2 2 / 3 
and k n = 2 n . Defining 

E n —u n u n 
H n = (2) 

for e = v e = there are two inviscid invariants E = J2E n and H = H n . The 
first is energy and the second, while not positive-definite, can be treated as a 
Hamiltonian with canonical variables u n and follows: 

du n jdt = ik n 5H/5u- n (3) 

The energy transfer between shells is e n = —k n _\A n . There is a trivial, unstable 
"Kolmogorov" solution of with u n = — i(2 1 / 3 e) 1 / 3 /c n ; 1 / 3 , H n = 0, and E n = 
22/9 e 2/3^,-2/3 correS p 0nc ii n g to a solution of an earlier cascade model || . 

The context for discussing this model along with the GOY model and 
Navier-Stokes is the extra invariant H. Despite the fact that this H is cu- 
bic and not quadratic, and that neither Euler nor GOY has a Hamiltonian 
of this form, due to the non-positive definite nature of H, it appears to have 
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some of the same qualitative effects upon the cascade that we speculate helicity 
is capable of for Navier-Stokes and GOY. The point is that while additional 
invariants can block the energy cascade, as enstrophy does in two dimensions 
and helicity does to some degree in Navier-Stokes j|, ||, due to the non-positive 
definite nature of the invariant there is an escape route where the cascade can 
find a way around this blockage. As will be demonstrated for the Kerr-Siggia 
model, this can take the form of pulses. From the tools used to demonstrate 
this for the Kerr-Siggia model, it will then be demonstrated that there is weak 
evidence for analogous phenomena in Navier-Stokes. 

In the original discussion 0, two classes of solutions besides the trivial "Kol- 
mogorov" solution were discussed. First, stationary solutions for a small num- 
ber of shells with no forcing or dissipation (e = v e — 0) and maximal H were 
discussed. Second, forced, dissipative solutions were discussed. Intermittency 
was found in the time dependent solutions and the effect of the extra invariant 
was noted, particularly as it affected the slope of the energy spectrum, which 
was < E n >~ k~ x l 2 rather than the Kolmogorov solution, but what effect the 
stationary solutions might exert upon the forced, dissipative time-dependent 
solutions was not considered. 

For the present calculations, e = 1 was chosen, which gives a characteristic 
timescale of t — 1. Using as initial conditions U\ = U2 = (1, 1), u n — 0, n > 3, 
N = 14, it took until t = 3.9 for the effects of initial transients to dissappear. 
Then statistics were taken until t = 6.8. Figure 1 shows the spectra of < E n > 
and < H n > for this period as well as k' 1 ! 2 curves, confirming the results of 
the original paper ||. Details will be discussed after the evidence for pulses is 
presented. 

Figure 2 shows E n and H n spectra for a series of moderately spaced times 
and the time development of E, H and dissipation for this time period. By 
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moderately spaced in time it is meant that the times shown are not so closely 
spaced so as to show continuous development, but are close enough to show 
a relationship between pulses of E n and H n and intermittent bursts in the 
dissipation. 

The primary event to focus upon is best illustrated in the H n spectra. In 
this sequence it starts as a pulse of positive H n centered on n — 4 at t = 6.11. 
It is associated with only one of several bumps in the energy spectrum at this 
time and is not associated with the spike in energy dissipation at t — 6.15. This 
spike in energy dissipation is assocated with one of the higher shell bumps in 
the energy spectrum and comes from a pulse at an earlier time of oppositely 
signed H n similar to the pulse to be described. 

Following the appearance of the positive peak of H n in n = 4 at t — 6.11, 
this peak breaks off from the lower shells and slowly propagates to larger shell 
numbers. The energy peak associated with it moves in tandem. Spectra of the 
transfer rates of E n and H n have also peaks that move with the pulse. When the 
effects of the highest shell, where the dissipation occurs, are felt, the pulse stalls 
at t = 6.29 before the energy in the pulse suddenly dissipates at t — 6.34. The 
stalling is the probable source of the bump in the time averaged energy spectra 
just before the dissipation regime. While this bump is on top of a spectrum 
less steep than Kolmogorov (k~ 1 / 2 rather than fc~ 2 / 3 ), it is qualitatively similar 



to a bump in the turbulent energy spectra for atmospheric observations ||10||, 
spectral closures^] and forced calculations of Navier-Stokes turbulence[p4|. For 
Navier-Stokes the bump is believed to be associated with a bottleneck effect p2 
where the decrease in the slope of kinetic energy spectrum in the dissipation 
regime blocks the free-flow of kinetic energy just at the boundary between the 
inertial and dissipation subranges. While this effect probably plays some role 
in the appearance of the bump in the Kerr-Siggia model, examination of figure 
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2 suggests a strong role for the stationary solutions associated with the second 
invariant. This comes from noticing that the cubic invariant is nearly maximal 
over the shells covered by the bump. 

While this pulse is dissipating at t — 6.34, the next major pulse of negative 
H n is beginning to move into shell 2 and positive H n for the major pulse fol- 
lowing that is developing in shell from the forcing. So a succession of E n and 
alternately signed H n pulses is suggested. Clearly this is a simplified picture as 
there are minor spikes in dissipation between the major spikes that are associ- 
ated with weak pulses with small H n of no particular sign. An example of such 
a weak pulse is the blip in H n at n = 8 for t = 6.35 and the rapidly moving 
E n at this time. To quantitatively demonstrate the alternation in sign of the 
strong pulses, figure 3 is a contour plot in shell and time separation of correla- 
tions between different shells and times of E n and H n . These plots are similar 



to contour plots of the energy transfer in forced Navier-Stokes calculations Hl2 



and also JT3[ and in meteorology are referred to as Hovmuller diagrams. These 
are: 



< {F n +A n ,t+At — F n+ A n )(F n<t — F n ) > (4) 

where F n is either E n or H n . Positive correlations are dark, negative are light. 
These plots are for n = 1, the second shell. The effect of a single pulse is the 
first region of increasing A n and At originating at (0,0). The propagation is 
linear after the first few shells. Starting at about At = 0.5 there is another 
strong dark region in the E n correlation and a strong light region in the H n 
correlation. This supports the qualitative picture coming from watching the 
time development that there are a succession of pulses of oppositely signed H n . 

The appearance of these pulses raises several questions. First, what modifi- 
cation of the stationary solutions can propagate as a unit? Second, what causes 
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the alternation in sign of H n of the pulses, is it the forcing or is it the nonlinear 
dynamics? We will not attempt to answer these questions. The point we do 
want to make is that there is some connection between the alternation in sign of 
the extra non-positive definite conserved quantity that seems to be associated 
with the appearance of pulses in the energy cascade and with intermittency in 
the model. In calculations where the extra invariant is suppressed, intermit- 
tency dissappears. These ideas are supported by noting that the mechanism 
with which the conserved quantities are pumped in the system and removed 



from the system can influence the scaling laws in the inertial subrange p5|, |26 



An extreme example is a calculation of the Kerr-Siggia model with a Newto- 



nian viscosity p5| where the extra conserved quantity is suppressed, there is a 
Kolmogorov spectrum and no intermittency. These are subtle questions that 
would require more accurate studies. 

How can the pulses be related to the spectra in figure 1? The < E n > 
spectrum in figure 1 goes as 2~ n / 2 . By dimensional arguments one might expect 
that the < H n > spectrum would obey 2~ 3n//4 , but this is not required since at 
any given time < H n > can have either sign. In fact, the < H n > spectrum is 
less steep than this and also seems to follow 2~ n / 2 . To understand this, imagine 
that each pulse is a coherent package of E n and H n traversing the spectrum, 
spending on average 2~ n / 2 time in each shell. Then the time averaged spectra 
< E n > amd < H n > will both have 2~ n / 2 spectra. This is similar to the 
argument that has been used to generate a -5/3 spectrum from fluctuations in 
a strained Burgers vortex^]. If At spent in each band goes as 2~ n / 2 as this 
suggests, then this would imply that the bands in the (A n , At) plots should 
approach zero slope as A n increases. There is some tendency in this direction 
for small A„ in figure 3, but for larger A n when the stalling noted at t = 6.29 
in figure 2 and dissipation effects are important, the bands are nearly linear. 
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Again, the time-averaged < H n > should not have any particular sign, as 
evidenced by shell 4 in figure 1, and their magnitude | < H n > | should 
decrease as the averaging time is lengthened. This has been verified by using 
different time intervals for the time averaging. 



3 The GOY model 

Given this discussion, let us now examine properties of the Kerr-Siggia model 
shared by the GOY model § 0. 

The GOY model has a very rich dynamical behaviour and it has been the 
object of many studies in recent years (see |7|, [15|, [16], [17], [L8|, |19| for some 



numerical and analytical results). It is the most popular shell-model for 3 di- 
mensional turbulence because of its intermittent properties are very close to 
the corresponding quantities in Navier-Stokes equations when the parameters 
of the nonlinear terms share some properties with the nonlinear term in the 
Navier-Stokes equations. In particular, for zero viscosity and no external forc- 
ing, when the system has the same conservation laws as a 3D flow: conservation 
of energy, of helicity and of volume in phase space. 
The dynamical equations follows: 

— u n = ik n (u n+l u n+2 + bu n+1 u n _ 1 + cu n _ 1 u n _ 2 ^j — vk 2 n u n + <5„ ino / (5) 

where v is the viscosity and / is a forcing acting on a large-scale-shell (for 
example, = 1) introduced to obtain a statistically stationary dynamical 
state. This model has interactions only between first and second-neighbor 
shells in the Fourier space. The two parameters b, c in the nonlinear terms are 
chosen such as to conserve energy, E = YJ n \u n \ 2 , for any choice of A. The most 
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general choice of parameters is: 

6 = c = -i^ (6) 

where e is the second free parameter in the model. 

The GOY model also has a second quadratic invariant beside energy: 

# = £x(e)X (£,A) kr- (7) 

n 

While energy conservation is forced by the choice (|^), the characteristics of the 
second invariant, H, change by changing the values of e and A. When e < 1 this 
second invariant is not positive-definite (x( e ) = — 1); while if e > 1 it is positive- 
definite (x( e ) — +!)• By remembering that the Navier-Stokes equations are 
characterized by having a second inviscid invariant that is positive-definite in 
2D (enstrophy) and non-positive definite in 3D (helicity), the value e = 1 can be 
identified as the border between a shell model for 2D turbulence (e > 1) and a 
shell model for 3D turbulence (e < 1). In the following, the problem of whether 



shell models like GOY-model are a good representation of 2D turbulence [pO 
is not addressed and only the range (0 < e < 1) where the dynamics should 
reproduce aspects of a 3D turbulent flow will be considered. 

By looking in detail at the structure of the second invariant, only when 
a(e, A) = 1 does it have physical dimensions coinciding with Navier-Stokes 
helicity J7|. This defines a line in the plane of free parameters where the 
inviscid conservation laws of the GOY model are very similar to the 3D Navier- 
Stokes equations. Because this invariant has the same physical dimensions as 
3D helicity and it is non-positive, we will denote it as the GOY -helicity in 
the following. In the last section, a modified version of the GOY model will 
be introduced with a second invariant having more correspondance with fluid 
helicity. 
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A necessary point before going on is that model (|5|) has two inviscid fixed 
points corresponding to the Kolmogorov scaling \u n \ ~ fc" 1 ' 3 (constant flux of 
energy, zero flux of helicity) and to a fluxless scaling \u n \ ~ k~( 1+a ^ 3 (constant 
flux of helicity, zero flux of energy) jn| . 

For this study our interest in this model comes from the presence in the (e, A) 
plane of a region where the static Kolmogorov-like fixed point is dynamically 



unstable |18| . The dynamics is fully chaotic and shows an intermittent cascade 
of energy toward small scales [|l^, [17]] with a complex (multifractal) structure of 
the attractor in the phase-space. This intermittency is quantified by measuring 
the scaling exponents ((p) for the structure functions in the inertial range: 

S p {k n ) =< \u n \ p >~ fc«P> (8) 

Only very recently [[?], [l£| has it been realized that the second quadratic in- 
variant plays a crucial role in the dynamics of the model. In it was found 
by varying the two free parameters (e, A) in the 3D-physically relevant region 
(0 < e < 1;A > 1) that along the line of constant helicity (a(e, A) = 1), the 
model has the same intermittent behaviour. That is, the set of ((p) depends 
only on the value of a, giving, for the first time, numerical evidence that the 
dynamics of the model is strongly dependent on the presence of the second 



inviscid-invariant. Furthermore, it has been shown [|19| | that by modifying the 
nonlinear term such as to destroy the presence of the second invariant (but 
still preserving the inviscid energy conservation) the intermittent corrections 
to K41 seem to weaken. 
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4 The transition to chaos 



It has been shown [|l8|j that by fixing A and varying the e parameters (and 
therefore, by changing a) the GOY model undergoes a transition to chaos fol- 
lowing a "Ruelle-Takens" scenario. In particular, there exists a critical value, 
e c such that for e < e c the Kolmogorov fixed point u n ~ k' 1 ' 3 is dynami- 
cally stable. Extending this analysis by changing the ratio between shells in 
the range 1 < A < 3 it is found that the "Ruelle-Takens" transition is quite 



general pip : there exists a line in the plane (e, A) which divides the region where 
a Kolmogorov-like fixed point is dynamically stable from a region where the 
dynamics is chaotic and intermittent. This qualitative trend of the transition 
line appears to beto be universal, even if the exact location can be slightly 
influenced by the forcing and by the value of viscosity. In the following, a very 
simple argument is presented based on the presence of the second non-positive 
invariant that predicts, with good accuracy, the existence of the transition and 
its location for any value in the plane (e, A). 

Consider the two inviscid quadratic invariants, the energy and the generalized- 
helicity, and their currents. First, for zero viscosity and zero forcing energy 
conservation gives: 

^^|^ra| Jn—1 Jni (9) 

where the energy current J n 

J n = 3f[-A n+1 -(l-e)A n ]. (10) 
is defined in terms of triple correlations: 

A n = kn-xUn-iUnUn+x. (11) 



The second conservation law for helicity takes the form: 

d 
dt 



— (— )"A; n |M n | 2 — L n _i — L n (12) 
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where the current of helicity, L n , from the nth shell to the n + 1th shell is: 



L n = (-) n k n S[A n -A n+1 }. (13) 

Let us suppose, for the moment, that there exists only one conserved quan- 
tity: energy. Then, very standard arguments |22| tell us that if viscosity is 
zero, the system tends to equipartition, corresponding in the GOY model to 
\u n \ 2 = const.. If one switches on viscous effects, and starts with an initial 
configuration with energy concentrated in the first shells, an energy cascade 
toward small scales develops. This energy cascade has been interpreted as the 
attempt of the system to reach new equipartition state |22|. This attempt at 



restoring equipartition is frustrated by viscous dissipation at small scales that 
continuously removes energy and prevents the small scales from reaching an 
equilibrium or quasi-equilibrium state. 

The Kolmogorov 1941 cascade, describing a smooth and constant transfer 
of energy from large scales to small scales is another way to rephrase this 
mechanism. But why does the flow not follow this picture of relaxing to a 
smooth and homogeneous transfer of energy and instead prefers to use a highly 
intermittent cascades consisting of bursts and blockages which are the origins 
of the intermittent corrections to the ((p) exponents? This is where we are 
proposing that the second inviscid quadratic quantity enters into the picture. 

It is well known that in 2D turbulence the presence of a second positive- 
definite quadratic invariant (enstrophy) does not allow the energy to cascade 
forward (toward small scales) [B3[]. This is a general result; either: the two 
conserved quantities must transfer in different directions in Fourier space or: 
only one can transfer to small scales. For example, in 2D turbulence it is widely 
believed that there exists a forward transfer of enstrophy and a backward 
transfer of energy (inverse cascade). 
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In contrast, the presence of a second non-positive definite quadratic invari- 
ant, like helicity in 3D turbulence, is only a minor constraint on the forward 
transfer of energy. Moreover, it is a constraint that, due to the non-positiveness, 
can have strong spatial and temporal fluctuations. If this picture is correct, 
intermittency in the 3D energy transfer could be the result of a competition 
between energy and helicity cascades. Temporal and/or spatial intermittency 
in the energy flux would be the result of switching between a net transfer of 
energy (possible due to cancellation effects in the helicity flux) and a depletion 
in the energy-transfer due to the presence of a non-zero helicity flux. 

How can these phenomenological ideas be checked in the GOY model? In 
the GOY model, a smooth and non-intermittent energy-transfer would corre- 
spond to dynamics near the Kolmogorov manifold u n ~ This implies 
that the energy flux (P) is almost constant in the inertial range and that the 
helicity flux (|12f) is almost vanishing. This Kolmogorov behaviour is obtained 
when the model has a static stable fixed point. It is natural, then, to ask if 
it is possible to understand the transition from the static behaviour to chaotic 
dynamics by invoking the second invariant. By plugging the Kolmogorov so- 
lutions into the expression for the generalized-helicity ([?]), we obtain: 



or negative determines whether H receives most of its important contribution 
from small or large scales, respectively. Therefore, when a > 2/3 the second 
invariant is concentrated at small scales and, as in 2D turbulence, prevents a 
smooth forward transfer of energy. This is reflected by strong intermittency 
and large deviations from Kolmogorov scaling. But, one can imagine that 
from time to time that it is still possible to transfer energy if some cancellation 




(14) 



n 



It is therefore clear that, whether the exponent 
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effects lead to an almost zero if -flux. On the other hand, when a < 2/3 energy 
transfers toward small scales without having any relevant change in H, i.e. the 
model relaxes in to a trivial Kolmogorov like fixed point. 

Figure 4 shows the numerical results J21J for the transition from a static 
Kolmogorov behaviour to chaotic dynamics by changing e and A. As predicted, 
the transition happens near the critical line defined by: 

«(e c ,A c ) = | A c = (1 - e c )- 3 / 2 (15) 

The systematic shift of 5% between the prediction (|i~5|) and the numerical 
results is probably due to viscosity as previously discussed |21|] We believe that 
this very simple result is a new important confirmation that the dynamics of 
the model is strongly influenced by the second invariant. 

Figure 5 shows decaying numerical experiments, that is zero forcing and 
non-zero viscosity with energy initially concentrated at large scales, that com- 
pare the dynamical behaviour of the model in two characteristic regions: case 
A [e = 0.1, A = 2], where there is smooth energy-transfer regime (on the right 
side of the critical curve in fig. 4) and case B [e = 0.5, A = 2], where the 
dynamics is chaotic and intermittent (left side of the curve in fig. 4). What is 
interesting is that for case A, when the second invariant does not introduce any 
constraint in the energy transfer (a < 2/3), the energy dissipation is a smooth 
function. This means that energy is transferred through the inertial subrange 
without any blocking and with a power-law behaviour in time. But, in case B 
(a > 2/3), the energy dissipation has a staircase shape in time indicating long 
periods when little energy reaches small scales (blocking) interrupted by short 
bursts of dissipation. This would be consistent with the suggested role for 
helicity where in addition to blocking the transfer, cancellation in the helicity 
transfer is associated with strong dissipation events. 
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Figure 6 shows the time evolution of the energy current J n and the helicity 
current L n through a shell in the inertial range where the model is chaotic. 
Obviously, the two fluxes are correlated, but the interesting fact is that when 
there is a burst of energy, the helicity flux has a sinusoidal shape, i.e. a net 
forward transfer of energy is only possible if the net averaged transfer of helicity 
is zero. Expanding upon the earlier proposal, this suggests blocking of energy 
transfer due to competition with helicity, interrupted by strong dissipation 
events made possible by brief, intermittent periods of large helicity fluctuations, 
but no net helicity flux. The dissipation events might then be associated with 
a strong dynamical coupling between modes with oppositely signed helicity, 
which would permit large helicity fluctuations, but no net helicity flux. 

5 A new shell model 

As described in the previous section, the structure of what we call helicity in 
the GOY model is only partially consistent with the helicity in the Navier- 
Stokes equations . Apart from the observation that it has the right dimensions 
and that it is not positive-definite, there is an asymmetry between odd and 
even shells that does not any counterpart in physical flows. One means of 
overcoming this problem is to introduce two dynamical variables in each shell, 
one transporting positive helicity w+ and the other transporting negative he- 
licity u~. The next step is choosing how to couple these terms. For this, 
we will use a complete decomposition of the three dimensional Navier-Stokes 
eqs. H into a basis where the two independent components of the velocity field 
at each wavenumber correspond to two pure helical waves. In such a basis 
there are 4 possible independent classes of triads interactions distinguished by 
the combination of helicity transported from each one of the three interacting 
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modes. 

Let us fix, for simplicity, the three modes q, p, k such that |k| < |p| < |q|, 
and call u Sk (k), u Sp (p), u Sq (q) the three interacting modes, where (sk,s p ,s q ) = 
(±1, ±1, ±1) refer to the sign of helicity in each mode. Then, it is simple to 
show that each triad can fall into one of the four following classes: 

1. (s fc , s p , s q ) = (+, +, +), or (-, -, -) 

2. (s k , s p , s q ) = (+, +, -), or (-, -, +) 

3. (s k , s p , s q ) = (+, -, -), or (-, +, +) 

4. (s k , s p , s q ) = (-, +, -), or (+, -, +) 

Following this decomposition one is led naturally to introduce 4 classes of 
shell models, each one corresponding to one of the four independent classes of 
triad interaction present in Navier-Stokes eqs. 

What is quite remarkable is that the original GOY model belongs to one of 
this classes (the fourth). To demonstrate this, we write the general equation 
for this class using positive helicity shells w+ and negative helicity shells u~. 

— u\ = ik n (u n+1 u^ +2 + bu n+1 u n _ 1 + cu n _ 1 u^_ 2 s j — vk n u\ + 6 n ,n f + 

jU~ = i k n (u+ +1 U- +2 + + CM+_l^n- 2 ) - + <Wo.T- (16) 

for which the conserved energy and helicity are given by: 

e = T,K\ 2 + K\ 2 (i7) 

n 

H = J2k n (\<\ 2 -K\ 2 ) (18) 

n 

exactly as in the Fourier- helicity decomposition of Navier-Stokes equation 0. 
By noticing that in the original GOY model shells n and n + 2 have the same 

17 



GOY-helicity, it can be seen that (|T6| ) is formed from two masked and mi- 
correlated versions of the original GOY model for the dynamical evolution of 
the variables (uf, , u£ , u^-i, u 2n) and u ti u^ n _ x , u^ n ). There- 

fore (^) has, by definition, the same behaviour as the previous model. From 
this, we think it is of primary importance to study in details the dynamical 
behaviour of the other three shell models (corresponding to the classes 1,2,3), 
allowing helicity to have all the dynamical interactions found in Navier-Stokes. 



Work in this direction is in progress and will be reported elsewhere [31 



6 Navier-Stokes helicity 

Up to this point only ideal shell models with extra non-positive definite invari- 
ants have been considered and how they might be extended to more closely 
resemble the helicity interactions in the Navier-Stokes equations. We would 
like to relate these ideas directly to the Navier-Stokes equations. As noted, 
earlier attempts at understanding the effects of helicity have emphasized its 
power to block the cascade 0, But, despite the blocking power of the extra 
invariant, the shell model calculations are indicating that the cascade can pro- 
ceed through interactions between shells where the sign of the extra invariant 
is opposite. In a full calculation, we would also want to see what the effects 
of helicity in physical space are. For example, there are low Reynolds number 
Navier-Stokes calculations of how vortex rings link and unlink and can generate 
and destroy helicity |29|, jSOfl . Because spectral properties were not analyzed and 
because of the low Reynolds numbers of these simulations, strong conclusions 
about the effects of helicity cannot be made from these calculations. But it can 
be said that even though the initial conditions contained large-scale helicity, 
small scale structures appeared and production of helicity from viscous effects 



18 



was not strongly blocked. 

Therefore, initial conditions that contain significant large-scale helicity, but 
show more clearly how a cascade is not blocked, would be desireable. Whether 
or not production of small-scale helicity by anisotropies plays a role, as has 
been suggested f28fl , will not be our objective. Simulations of isotropic, homo- 
geneous turbulence in a periodic box are traditionally initialized with a given 
spectrum, but the phases of individual wavenumbers is completely random. To 
test the effects of helicity we propose constraining the phases of the velocity 
components such that the helicity of the Fourier modes is not all of one sign, 
unlike the investigations of the blocking power of helicity ||. Since shell models 
are indicating strong effects from extra invariant fluctuations even when the 
average value is zero, this suggests initializing a Navier-Stokes calculation with 
net zero helicity. 

Several tests of this type have been done, all with one qualitatively similar 
feature. Helicity when it first appears in the spectra pops up in two neighboring 
bands of opposite sign, then the bands separate. Other than this qualitatively 
feature, the number of tests is not yet sufficiently large to make definitive 
statements. Two cases are shown in figure 7. Both cases are 64 3 simulations 
where only a small number of modes in wavenumber band 4 were initialized. 
Each has the common feature noted, but a rich variety of additional features 
as well. 

In case A, each mode is initialized with maximal helicity, but otherwise 
the phases were chosen randomly (by hand) and the net helicity was zero. 
Very quickly the helicity picture changes. From a zero helicity spectrum, soon 
helicity of opposite sign appears in shells on opposite sides of the initial energy 
shell. For a short period, time sequences show that the helicity peak at higher 
wavenumber moves to small scales until it dissipates, leaving net helicity of the 
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sign of the large scale peak. It is during this phase that relative dissipation 
rate (dissipation/energy) is largest. Therefore, through dissipation of helicity 
at small scales, large-scale helicity is generated. Once only the large-scale 
helicity peak is left, the blocking action of the helicity at large-scales becomes 
important and the dissipation rate is suppressed. 

In case D, the helicity of each mode is chosen to be zero by using free- 
slip boundary conditions in along central planes in the box. All the modes 
except one use the same free-slip symmetry plane, with the exception imposed 
to break this symmetry. From these initial conditions helicity does not initial 
grow around the initial k — 4 wavenumber band, but around the resonance 
band at k — 8. Note that once again helicity first appears in neighboring bands 
with opposite sign. Then time sequences show that the bands move towards 
opposite ends of the spectrum. Once again, dissipation is largest when the high 
wavenumber helicity band moves into the dissipation band and is annihilated, 
then decreases when large-scale helicity of only one sign remains. 

To make this more quantitative, more calculations need to be done and there 
needs to be analysis of helicity spectra and transfer properties. But what of a 
relationship to the shell models? Investigations of shell models with many more 
than 2 variables per shell seem to invariably lead to reductions in intermittency. 
Fully developed turbulence when viewed as a shell model has an infinite number 
of degrees of freedom and should in this sense not be intermittent. But it is, and 
this is understood as being due to coherent structures in physical space, which 
in Fourier space implies strong phase correlations. The phase correlations 
therefore prune the number of paths the cascade can take, returning us to 
simple models with few paths and strong intermittency, such as those presented 
here. Therefore, to get meaningful comparisons between 3D direct calculations 
and shell models, there must be some means of identifying the paths the cascade 
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will follow and calculating statistics along these paths. Given the difficulty of 
attaining this, let us make some other suggestions. 

First, there needs to be further work on bi- dimensional correlations of 
wavenumber and time, similar to figure 3 here and in earlier analysis [|12|, |I~3" 



The question with direct calculations is what quantities to use. It has been 



found |L2| that energy transfer spectra have a strong signature. Helicity spec- 
tra and helicity transfer spectra need to be analyzed in the same manner. For 
simulations with a small number of initial modes, such as the examples just 
given, at least for short times statistics for modes formed by the initial interac- 
tions and their daughters could be studied. It is our hope that analysis of new 
Navier-Stokes simulations of this type, coupled with any new understanding fo 
the role of helicity coming from shell models, will provide new insight into the 
nature of the intermittent cascade of energy to small scales in turbulent flows. 
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FIGURE CAPTIONS 



1. Time-averaged spectra for the kinetic energy (solid line) and Hamiltonian 
(dashes) for the Kerr-Siggia model. Note the "bump" in the energy 
spectra near shell 9. The gap in the Hamiltonian spectra at shell 4 is 
where it is negative. 2~ n l 2 curves are drawn for comparison. 

2. Spectra at individual times for the kinetic energy and Hamiltonian. In 
order, the times correspond to 1st time: — , 2nd time: -, and 3rd 
time: ■ • •. Also given is the time evolution of the total kinetic energy — , 
Hamiltonian — and energy dissipation • • • for the period covered by the 
spectra. The X in the energy curves indicates roughly the times of the 
spectra to the left. 

3. Hovmiiller diagrams in shell and time separation of kinetic energy and 
Hamiltonian fluctuations at different shell numbers and times using (f|). 
Negative fluctuations correlations are light, positive are dark. Note the 
strongest black for energy is not at (0,0), indicating that the strongest 
fluctuations in energy are in the higher shells. The sign of H fluctuations 
changes between pulses (light upper ridge), but does not for energy (dark 
upper ridge). 

4. Comparision in the (e, A) plane of the, numerically estimated, transition 
(circles) |21j and the theoretical prediction (solid line). 



5. Log-log plot of the total energy as a function of time in a pure decaying 
simulation. Case A (dotted line) corresponds to a smooth Kolmogorov- 
like transfer of energy, while case B (solid line) corresponds to a chaotic 
intermittent energy-transfer. 
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6. Fluxes of energy and helicity during a burst trough a shell in the inertial 
range. Notice the oscillatory behavior of the helicity flux triggering the 
energy transfer. 

7. Navier-Stokes three-dimensional banded wavespectra of kinetic energy 
and helicity at individual times. At t — only energy is shown for each 
case. For case A, helicity at later time (t = 0.75) is multiplied by 10. 
For case D, helicity at earlier time (t = 0.0625) is multiplied by 10. 
Triangles indicate helicity maxima at the two times for each case and 
crosses indicate helicity minima. 



26 



